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We evaluate the neutrino nucleon charged current cross section at next-to-leading order in quan- 
tum chromodynamic corrections in the variable flavor number scheme and the fixed flavor number 
scheme, taking into account quark masses. The number scheme dependence is largest at the highest 
energies considered here, 10 12 GeV, where the cross sections differ by approximately 13%. We illus- 
trate the numerical implications of the inconsistent application of the fixed flavor number scheme. 



I. INTRODUCTION 

Inelastic scattering experiments with lcptons interact- 
ing with proton and neutron targets and hadron-hadron 
collider results have provided a picture of the structure 
of the nucleon over a wide range of momentum transfers. 
These data combined have led to a parton model picture 
with parton distribution functions (PDFs) extracted by 
a number of groups [ll-@ . With large underground neu- 
trino telescopes designed to detect neutrinos from astro- 
physical sources |7Hl0f . the PDF inputs to the neutrino- 
nucleon cross section are essential ingredients to uncover 
features of the astrophysical sources. 

Heavy quark masses and their roles in the theory of 
the structure functions and the extraction of parton dis- 
tribution functions at next-to-leading order have been ex- 
plored extensively [TT| - [25| . As neutrino telescope analy- 
ses become more refined, it is useful to consider the quark 
mass effects in the evaluation of neutrino-nucleon cross 
sections. Furthermore, in considering ultrahigh ener- 
gies, one can explore the implications of potentially large 
\og{Q 2 /tuq) corrections from non-zero quark masses. 

We use the neutrino-nucleon charged current scat- 
tering example to discuss two different theoretical ap- 
proaches to including heavy quark contributions, in- 
cluding next-to-leading order quantum chromodynamic 
(NLO QCD) corrections. As the incident neutrino en- 
ergy increases, the average momentum transfer Q also 
increases, so heavy quarks become effectively "light" fla- 
vors. We look at the fixed flavor number scheme (FFNS) 
where the number of light quark flavors is fixed, regard- 
less of the scale of Q 2 . The Gluck, Reya and Vogt 
(GRV) [26| 3-flavor PDFs are useful for this scheme. 
Gluck, Jimenez-Delgado and Reya (GJR) |4| have up- 
dated these 3-flavor PDFs. The variable flavor number 
scheme (VFNS) allows for the introduction of charm and 
bottom quarks as constituents of the nucleon as Q 2 in- 
creases. We use the CTEQ6.6M PDFs Q and a version 
of the GJR PDFs [E[ which are applicable in the VFNS. 
The GJR variable flavor PDF set is generated radiatively 
from the 3-flavor set at a factorization scale set by the 
heavy quark mass tuq . The GJR variable flavor PDF set 
is not fit to data beyond the initial three-flavor fit. 

The modified minimal subtraction scheme for NLO 
corrections (MS) with massless quarks is a variable flavor 
scheme which neglects the quark mass except effectively 
with a step-function threshold factor. A version of the 



Avizas, Collins, Olness and Tung prescription (ACOT) 
[ilj . fl9| is the approach to include heavy flavor [2l| in 
the variable flavor number scheme that we use here, al- 
though there are other approaches to incorpora ting the 
quark mass effects in the generalized mass VFNS |22h25I | . 

In Section II, we review the formalism for the neutrino- 
nucleon charged current cross section at leading order 
and next-to- leading order in QCD, and the ACOT for- 
malism for including quark mass effects is described. At 
ultrahigh energies, the parton distribution functions are 
probed at very small momentum fractions, in some cases 
beyond the PDF fits. We discuss our extrapolation to 
smaller momentum fractions in Section II. 

Our results for the the NLO QCD corrected neutrino- 
nucleon charged current cross sections are shown in Sec. 
III. The focus in this paper is on the high energy regime, 
however, we also show the 100 GeV to 10 4 GeV energy 
regime, where one can see effects of flavor scheme, pre- 
scription choice and PDF choice. The NLO QCD cor- 
rection is on the order of 3% in this energy regime. At 
ultrahigh energies, quark mass effects (with the exception 
of the top quark) are negligible so the massless MS for- 
malism represents the VFNS. The VFNS accounts for a 
resummation of logarithms ln(Q 2 /mg) which the FFNS 
neglects. 

As we see numerically in Sec. Ill, the FFNS with only 
three light quarks yields an enhancemnt of the VFNS 
neutrino-nucleon charged current cross section by more 
than 10%. Already by neutrino energies of ~ 10 7 GeV, 
discrepancies appear in the cross section calculated with 
and without the resummation of these logarithms. We 
summarize our results in Sec. IV, where we also show 
the numerical implications of a mismatch in scheme at 
ultrahigh energies p7| . 



II. NEUTRINO NUCLEON CHARGED 
CURRENT SCATTERING 

A. Leading-order cross section 

In neutrino-nucleon charged current scattering, one is 
interested in the inclusive process. We evaluate the cross 
section for neutrino scattering with isoscalar nucleons 
N = (n + p)/2. For definiteness, we consider muon neu- 
trino scattering. For the energies of interest here, we 
neglect the muon mass. The effect of the tau mass in 
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v T N scattering is on the order of 5% at E u = 10 3 GeV 
[28l | . so the "low energy" results reported here are not 
applicable to tau neutrino charged current scattering. 
With the momentum assignments, 
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In the massless quark and massless target limit, x is the 
momentum fraction of the nucleon carried by the struck 
parton. 

The neutrino cross section, in terms of the structure 
functions Fx, F% and F3, is 
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where Mw is the mass of the W boson and Mm is the 
mass of the target nucleon. The upper sign in Eq. ([5]) is 
for neutrinos, the lower sign for anti-neutrino scattering. 

At leading order, the quark mass and nucleon mass 
introduce two types of corrections which appear in the 
structure functions. In the massless limit, the light cone 
momentum fraction in the PDFs is the momentum frac- 
tion x. Including mass effects, the light cone momentum 
fraction changes to 
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yja 2 + b 2 + c 2 - 2(ab + bc + ac) (8) 



in terms of initial quark mass m\ and final quark mass 
m,2- The light cone momentum fraction £ goes into the 
quark distribution function evaluation. We have kept 
the nucleon mass correction in £ (through rj) and in the 
differential cross section. Below E v — 100 GeV, more 
target mass effects should be included (29j . 

A variant to this approach is to include the kinematic 
suppression associated with quark masses is through the 
variable x defined by [2(| 



x = m 1 + 
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For neutral current scattering of a neutrino with a c, this 
accounts for the fact that an associated c is also a compo- 
nent of the nucleon, leading to \c = Tj(l + 4m 2 /Q 2 ). Eq. 



([9|) extends this to charged current interactions where fi- 
nal quark and the remaining sea component quark masses 
are different. We discuss the numerical difference be- 
tween using £ and x i n t ne NLO cross section below. 

The second contribution comes in the mass corrections 
to the structure functions in terms of the quark distribu- 
tions. At leading order, neglecting corrections propor- 
tional to 1 + 4x 2 M 2 1 Q 2 , the mass corrections are [l8j 
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for an initial quark qi which converts to quark qj or the 
corresponding antiquark initiated process. The quantity 
Vij is the element of the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix. We use the central values quoted by the 
Particle Data Group [3(| ■ For neutrino scattering, i — 
d, s, b, j = u, c in the variable flavor number scheme. 
In the fixed flavor number scheme, i = d, s, j — u for 
the PDFs, however all flavors are included in the final 
state sum for the CKM matrix elements. The quantity 
\x is the factorization scale, which we set to /1 = Q in our 
numerical evaluations. 



B. Next-to- leading order corrections 



At next-to- leading order in QCD, the structure func- 
tions Fi have corrections which account for perturba- 
tive loop corrections and splitting of quarks, antiquarks 
and gluons in the nucleon. Graphically, the loop correc- 
tions and the quark splitting corrections are shown in 
Fig. [TJ In addition to the quark initiated contributions 
shown in this figure, there are also antiquark initiated 
diagrams. The gluon contribution to the NLO neutrino- 
nucleon cross section comes from the graphs in Fig. [2] 

The final evaluation of the structure functions and ul- 
timately the cross section requires a subtraction for fac- 
torization [3l| . The ACOT scheme uses quark masses to 
regulate collinear divergences. The proof of factorization 
by Collins in Ref. 31] is valid for any scale Q relative 
to mi. The ACOT prescription, for example at NLO 
for the gluon splitting in Fig. [5] with W*G — > be, has 
subtraction terms proportional to the gluon splitting co- 
efficient Pc-^q times In /j, 2 /m 2 and Pc-> q times In /j, 2 /m 2 . 
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FIG. 1: The NLO loop correction to the structure functions 
comes from the interference of the two upper diagrams dia- 
grams. An additional correction comes from the two lower 
diagrams. 





FIG. 2: Gluon splitting produces a quark-antiquark pair. For 
the fixed flavor number scheme with Nf = 3, this process is 
the only way to get charm quark and b quark contributions 
to the cross section. 



Generically, for example for the structure function Fi, 
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in terms of convolution integrals denoted by the symbol 

For massless quarks in the modified minimal subtrac- 
tion scheme MS, the logs are replaced by 1/e where di- 
mensional regularization (d — 4 — e) is used to regulate 
the infrared and collinear divergences. Even using small 
quark masses in eq. (11) for the gluon terms (g — > q%qj) 
and the massless MS for the quark terms (the S-ACOT 
scheme), the numerical results using massless MS are re- 
produced. With the ACOT and S-ACOT approaches 



where the charm and bottom quark masses are explicit, 
one can use the same formalism for a range of energy 
scales as the role of the quark changes from "heavy" to 
"light" 32]. Of course, the up, down and strange quarks 
are always "light" quarks in our evaluation. 

We use the ACOT(x) (or S-ACOT(x), as labeled) pre- 
scription [2(| for the inclusion of all quark masses. In the 
ACOT prescription for VFNS evaluations, the minimum 
z depends on whether or not the PDF in the convolution 
is for a gluon or quark/antiquark. The limits are 
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where £ and r\ are defined in eqs. (6) and (7). In the sub- 
traction terms, the massless splitting functions are used, 
however masses are kept in the coefficient functions oj. 
Eq. (|11[) shows general terms for the splitting functions 
and coefficient functions, but we note that at leading or- 
der, e.g., Pi_n = P]^] m the massless quark limit. The 
ACOT(x) prescription replaces £ — >• x m both the inte- 
gration limits and in the PDF of the leading order term. 

The large logarithms associated with the quark mass 
terms can cause numerical issues at high energies. We 
have made a series of approximations in our evaluation 
of the variable flavor number scheme to avoid these nu- 
merical problems. We keep all the masses at low energies, 
and at the highest energies, we evaluate the cross section 
with zero quark masses. 

The S-ACOT(x) prescription is used for neutrino en- 
ergies above E v — 10 6 GeV. As noted above, the S- 
ACOT(x) prescription uses the massless MS scheme for 
the quark initiated terms. The subtraction term for the 
gluon splitting to quarks also has the quark masses set 
to zero, however, the quark masses are retained in the 
gluon splitting process [32]. For E v > 10 9 GeV, we use 
the massless MS scheme, that is, the zero-mass version 
of the variable flavor number scheme. These approxima- 
tions are reliable for a wider range of energies than used 
here 

The ACOT scheme and its variants are only one in 
a class of approaches [22T - [24l | to including quark mass 
corrections in a general mass, variable flavor number 
context. One alternative is the Thorne- Roberts (TR) 
method which uses Q 2 = ttiq as a transition point for 
matching fixed flavor structure functions below the heavy 
quark mass (toq) to variable flavor structure functions 
above the transition point [22|. Another alternative is the 
FONLL method applied to DIS (H)]. The various meth- 
ods differ by terms which are subleading, but which nev- 
ertheless may have phenomenological implications. Ref. 
24] considers a rescaling alternative to x- A detailed 
comparison of the heavy flavor contributions in differ- 
ent generalized mass variable flavor number schemes to 
heavy flavor contributions to deep-inelastic electropro- 
duction structure functions Fi and Fl , with and without 
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X-scaling, appears in Sec. 22 of Ref. f25j. 

In Sec. Ill below, we focus on the ACOT(x) scheme 
and its variants, and a comparison of the high energy 
neutrino cross section evaluated with this specific gen- 
eralized mass variable flavor number scheme to a flavor 
number scheme fixed at low energies. Our primary focus 
is on how the cross sections differ at high energies, al- 
though we show results for energies as low as E v = 100 
GeV. We comment below on the magnitude of the vari- 
ations in the ACOT schemes (ACOT, S-ACOT and the 
X variants) in the 100 GeV range as a rough order of 
magnitude of the more general scheme dependence. A 
full comparison of the different generalized mass, vari- 
able flavor number schemes is beyond the scope of this 
paper. 

The comparison of a VFNS with the FFNS at high 
energies is useful to determine the numerical effects of 
the VFNS's resummation of \n(Q 2 /nig). In the fixed fla- 
vor number scheme with flavor number equal to three, 
heavy flavors and light flavors have a separate treatment 
regardless of the relation between Q 2 and characteristic 
m 2 q . The Gluck, Jimenez-Delgado and ReyaQ fixed fla- 
vor number scheme PDFs with three light flavors (u, d, s) 
carried in the evolution of the parton distribution func- 
tions are used here. We denote the FFNS version of the 
GJR PDFs by GJRF. Heavy quark contributions (c, b) 
appear only through explicit contributions from gluon 
splitting. In eq. (|11[) , in the fixed flavor number scheme 
with Nf = 3, the heavy quark contributions come in only 

through a/ 1 ' .-. There are no subtraction terms corre- 
sponding to heavy quarks since there are no heavy quark 
PDFs. 

For variable flavor number scheme evaluations, we use 
the GJR VFNS version Q in which the heavy quark con- 
stituents are radiatively generated from the 3 flavor fits to 
data, denoted GJRV. We also use the CTEQ6 [l| version 
incorporating heavy quark effects, CTEQ6.6M PDFs 
which updates the CTEQ6HQ version incorporatingthe 
quark mass effects through the ACOT prescription 

C. Extrapolation to small x 

As we are considering neutrino energies up to the high- 
est energy cosmic rays, E ~ 10 12 GeV, extrapolations of 
the PDFs beyond the range of experiments is required. 
The PDFs are available numerically for a range of Q 2 and 
for x m i n < x < 1. The characteristic Q 2 for high energy 
neutrino-nucleon scattering is Q 2 ~ since the prop- 
agator suppression dominates the evolution of the PDF 
with increasing Q 2 (34[. This is well within the range for 
both the CTEQ and GJR PDFs. The minimum value of 
x is more constraining at ultrahigh energies. Using 

xy(2ME v ) = Q 2 

and (y) ~ 0.2, x is required below x^® 6 = 10~ 8 for 
CTEQ6.6M and below x^ = 10~ 9 for GJR PDFs. 



We extrapolate the PDFs below x m ; n using a power 
law extrapolation [3^, [36| , where 

Xq(x : Q ) ^min^?(^min; Q )(*^/*^-min) 9 

xg(x,Q 2 ) = x min g(x min ,Q 2 )(x/x min y 3 . (14) 

The antiquark distribution equals the sea quark distribu- 
tion, and at low x, the valence contribution is negligible. 
For the CTEQ6.6M PDFs, we have used a log(Q) de- 
pendent form for Ai, while for the GJR sets, we use a 
constant A^. The values used for A; are shown in Table 
UJ where for the CTEQ6.6M set, we show the value for 
Q = M w . 
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TABLE I: The parameters that appear in the small x extrap- 
olations of eq. iffl]). For the CTEQ6.6M set, we use a log(Q) 
dependent form. Here we show the value at Q = Mw- 

Extrapolations using functional forms other than 
power laws have been suggested by a number of authors 
[37l f38j| . The typical range of predictions for E v = 10 12 
GeV is on the order of a factor of 0.5 — 2 times the cross 
sections reported here. 

III. RESULTS 

In our evaluation of the quark mass effect on the 
vN cross section at NLO, we restrict our attention to 
the CTEQ6.6M and GJR PDFs. We set the heavy 
quark masses to m c = 1.3 GeV, and mj= 4.5 GeV for 
CTEQ6.6M, and to m c = 1.3 GeV, and m b = 4.2 GeV for 
GJR as indicated in each PDF. As noted above, we use 
the ACOT(x) scheme for our evaluations in the variable 
flavor number scheme, with transitions to S-ACOT(x) 
and MS as the neutrino energy increases. 

In Fig. [21 we show the MS neutrino-nucleon charged 
current cross section as a function of the incident neutrino 
energy. This sets the scale of the cross section. The 
CTEQ6.6M PDFs give a cross section, using the massless 
MS scheme, that are within 2% of the standard CTEQ6 
result, using the MS scheme at E v = 100 GeV, and less 
than 0.5% different at E v = 10 10 GeV. In subsequent 
figures, we show only ratios. 

To get an idea of the relative importance of the fla- 
vor components, in Fig. [4] we show the ratio of the 
cross section for each sub-process in the neutrino nucleon 
charged-current cross section. At NLO, it is not possi- 
ble to separate, e.g., the d and u contributions since the 
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FIG. 3: The vN cross section for the charged current pro- 
cess as a function of the incident neutrino energy. The cross 
section is evaluated using the massless MS scheme with the 
CTEQ6.6M PDFs. 
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FIG. 4: The ratio of the separate flavor contributions to the 
ncutrino-nucleon charged current cross section, evaluated us- 
ing the CTEQ6.6M PDFs in the VFNS. 



gluon splitting diagrams (Fig. [2]) are added at the am- 
plitude level. The figure shows that the du contribution 
dominates until the charm mass corrections and valence 
contributions are neglible. Then, the sc contribution is 
nearly equal to the du contribution. At high energies, the 
su and dc terms are also nearly equal. At low energies, 
the valence d component more than compensates for the 
mass suppressed charm quark production. Contributions 
involving the b quark are at most at the level of 0.1%. 

Fig. [5] shows KNLO-factor, which is the ratio of NLO 
cross section to LO cross section, for the incident neu- 
trino energy between 10 2 GeV and 10 4 GeV. In each 
ratio, the LO cross section is evaluated using the same 
PDF as the NLO cross section, namely the NLO PDF 
set, to exhibit the size of the partonic cross section cor- 
rection. We also use x as the scaling variable for all but 
the massless MS result. A comparision of the different 
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FIG. 5: KNLO-factor, the ratio of the (tnlo to olo (with 
quark masses) for the VFNS and FFNS. The LO cross section 
is evaluated using the same PDFs as the NLO cross section 
and the scaling variable x ls used throughout. For VFNS, all 
quark masses are kept while for the MS curve, quark masses 
are set to zero. 



schemes and PDFs shows that at 100 GeV, they differ by 
as much as ~ 4%. The K NLO -factor of the GJRV VFNS 
with all masses included are very close to the massless 
MS results for E v > 10 3 GeV, where the quark masses 
have little impact. Below this energy, quark mass effects 
suppress some of the QCD corrections. The CTEQ6.6M 
results are intermediate between the massless MS GJRV 
and GJRV massive VFNS results below ~ 1 TeV. 

The KNLO-factor for the FFNS is only a little higher 
than the VFNS at E v = 10 2 GeV, but the difference in 
the KNLO-factor increases to about 1.5% relative to the 
GJRV KNLO-factor at 10 4 GeV. Since the KNLO-factor 
is a ratio, Fig. [5] does not illustrate the fact that the 3 
flavor NLO charged current cross section using GJRF is 
about 1% lower than the GJRV cross section. The GJRF 
"LO" cross section is lower than the GJRV "LO" cross 
section by ~ 2.5% at E v = 10 4 GeV (using the NLO 
PDFs). 

Fig. [5] also shows the difference in Kmlo between the 
ACOT(x) and the S-ACOT(v) prescriptions using the 
CTEQ6.6M PDFs at E v = 100 GeV, amounting to more 
than 3%. A similar positive offset of 3% is seen when 
one compares the cross section at NLO using the ACOT 
or S-ACOT scaling variable £ rather than the variable x- 
At 100 GeV, this would indicate that the specific vari- 
able flavor scheme for incorporating the quark mass and 
the scaling variable are more important than the fixed 
flavor/variable flavor choice. At higher energies (above 
1 TeV), ACOT and S-ACOT NLO cross sections differ 
by less than 1%, with the difference between using £ and 
X for the quark terms less than 1.5%. At E u = 10 4 
GeV, the NLO cross sections for the CTEQ6.6M PDF 
are essentially identical for ACOT, S-ACOT, ACOT(x) 
and S-ACOT(x). The differences between generalized 
mass, variable flavor number schemes and scaling vari- 
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able choices continue to be areas for further work [2£ 
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FIG. 6: KNLO-factor: the ratio of the neutrino-nucleon 
charged current <jnlo to oxo for E v = 10 2 — 10 12 GeV. The 
LO cross section is evaluated using the same PDFs as the NLO 
cross section. The VFNS cross section is evaluated using the 
ACOT(x) or S-ACOT(x) scheme for E„ < 10 9 GeV. 

In Fig. [5J we show the KNLO-factor for the full energy 
range for the GJR PDFs. As shown in Fig. [SJ the KNLO- 
factor for the VFNS and MS is essentially identical above 
E v ~ 10 3 GeV. Given the broad energy range of concur- 
rence with VFNS with masses and the massless MS, we 
use the massless MS cross section for the VFNS above 
E v = 10 9 GeV to avoid numerical errors associated with 
the subtraction terms. This "patching" is used for the 
remaining figures. 

The KNLO-factor for the FFNS starts to deviate from 
the VFNS K N LO-iactor at about E v ~ 10 4 GeV, with sig- 
nificant deviations by E v = 10 6 — 10 7 GeV where, in Fig. 
21 the charm quark contribution is effectively "massless." 
For evaluations at the level of less than 5% error in the 
neutrino nucleon cross section, this is the energy range 
above which the VFNS should be used. This is shown 
graphically in Fig. [7J 

In Fig. Efa) we show the ratio of VFNS NLO cross 
section to FFNS result for different PDFs. The GJRV 
set is used for the VFNS and the GJRF set for the FFNS 
result for the dashed GJR curve. The CTEQ6.6M set 
includes 5 quark consituents, so the "FFNS" in this figure 
for CTEQ6.6M simply omits contributions for c and b 
quarks and antiquarks, even though they are consituents 
of the nucleon in this set. The GJR sets are a better pair 
of PDFs to compare, since the two sets are designed to 
accommodate different flavor numbers. 

While the ratio of VFNS to FFNS cross sections in 
Fig. [TJa) are stable for energies below ~ f 6 GeV for the 
GJR PDFs, the ratio decreases as the incident energy is 
increased higher, and there is about f 3% discrepancy at 
I0 12 GeV. For CTEQ6.6M PDFs, they have almost the 
same value up to 10 5 GeV, and at higher energies their 
ratio varies depending on the energy. Their maximum 
difference is at highest energy, which is about 3%. 
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FIG. 7: (a) The ratio of the NLO charged current ctvfns 
to a FFNS for CTEQ6.6M and GJR PDFs. (b) Comparison 
of GJR and CTEQ6.6M PDFs for the NLO charged current 
ctvfns and uffns- 



We also compared the cross sections of the VFNS and 
FFNS for GJR and CTEQ6.6M PDFs. As shown in Fig. 
Ulb), the FFNS evaluations using the CTEQ6.6M and 
GJRF PDFs have very close values up to E v ~ I0 6 GeV. 
At higher energies, however, the FFNS cross section for 
GJR PDFs grows as the energy increases, and it makes 
difference about 15% with the result for CTEQ6.6M. The 
more reliable VFNS results for GJRV and CTEQ6.6M 
differ by about 5% at E v ~ 10 6 GeV, and differ by less 
than 3% at 10 12 GeV. We note that direct measurements 
of deep-inelastic scattering have been done only to an 
equivalent neutrino energy of E v ~ 5 x 10 4 GeV [39| . Up 
to this energy, the FFNS and VFNS ratios are essentially 
unity. 



IV. CONCLUSIONS 

A concurrence between number scheme, PDF set and 
related subtraction terms has been emphasized in the 
literature, e.g., in Ref. [25|, [33|. Mismatches in appli- 
cation can lead to errors on the order of 20%. To ex- 
hibit one such mismatch, we show the KNLO-factor for 
the CTEQ6.6M PDFs in Fig. |U The dot-dashed line, la- 
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beled "mixed scheme," shows the ratio of the NLO cross 
section to the LO cross section, where the NLO eval- 
uation has only d, s and u light quarks but with the 
full 5-flavor MS gluon subtraction correction. A similar 
KNLO-factor appears in Ref. [271 ]. in which the 3- flavor 
GRV PDFs are used. 

At energies below E v ~ 10 6 GeV, the ratio of the 
VFNS and FFNS cross sections using the GJR PDFs is 
~ 1, but at higher energies, the ratio drops. The almost 
13% discrepancy between the two cross sections can be 
attributed, at least in part, to the summation of large 
log(Q 2 /m|) in the VFNS PDFs. This is a quantitative 
example of the statement that a three flavor calculation 
of structure functions overestimates the "true" structure 
functions when more flavors are active [33] ]. 

The GJRV PDFs are not fit to data, but instead gen- 
erated radiatively from the 3-flavor fit. The CTEQ6.6M 
set is fit including mass effects, so they should be con- 
sidered more reliable at high energies, especially where 
the charm quark contribution is more important. Even 
so, the discrepancy between the GJRV and CTEQ6.6M 
VFNS evaluations agree well at the highest energies. 

Cooper-Sarkar and Sarkar (CSS) have evaluated the 
NLO neutrino nucleon cross sections using an indepen- 
dent fit to the data at NLO [HJ. Our vN charged 
current cross sections using CTEQ6.6M are bigger than 
CSS's cross sections by about 8-18% for s = 2M^E V = 
10 8 — 10 12 GeV 2 , generally within their estimates of PDF 
uncertainty. The GJRV cross section at s = 10 12 GeV 2 
differs from CSS by about 20%. 

Uncertainties at the level of a few or a few tens of per- 
cent at E v = 10 12 GeV rely on perturbative QCD and 
DGLAP evolution being applicable to very small x values 
for Q 2 ~ M w . At lower values of Q, gluon recombina- 
tion and saturation effects are important [37j ]. however, 
at Q = Mw, it is not clear that saturation should be 
important for the total cross section [4l| . As noted in 
the introduction, there are a range of predictions that 
do not rely on DGLAP evolved PDFs [38[. On a short 
time scale, one looks forward to further information as 
the LHC analyses yield PDFs from data in new ranges of 
x for Q ~ Mw, as a start to the experimental probe of 
PDFs and structure functions required for the ultrahigh 
energy neutrino cross sections. 
Note added: 

Since our submission of t his p aper, Gluck, Jimenez- 
Delgado and Reya, in Ref. [42| have emphasized that 
the K-factor is traditionally defined as the ratio between 



NLO and LO (with LO partonic cross sections and LO 
PDFs). To avoid confusion, we have relabeled our ratio 
from K in the original version of this paper to Knlo- 
We confirm the results of Ref. [42j that the K-factor 
as traditionally defined docs decline to about ~ 0.6 for 
the VFNS (GJRV) as compared to ~ 0.8 for the FFNS 
(GJRF). The authors of Ref. [42[ also comment on the 
importance of the b — t contribution at ultrahigh energies. 
Our conclusions about the ratio of VFNS/FFNS at NLO 
with the GJR PDFs do not change with the inclusion 
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charged current onlo to glo for E v = 10 — 10 GeV. The 
CTEQ6.6M PDFs are used. For the FFNS NLO results, the c 
and b PDFs are set to zero, but all five flavors of CTEQ6.6M 
are used for the LO cross section used for all the curves in 
the figure. The "mixed scheme" ratio uses 3 flavors of quark 
PDFs but makes a subtraction for 5 flavors from the gluon 
fusion term. 



of the top quark contribution. With the b — t term, the 
CTEQ6.6M charged current cross section is larger than 
the CSS cross sections by 13-31% for s = 10 8 - 10 12 
GeV 2 . 
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